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Abstract 

Ferroelectricity has been found to occur in several insulating systems, such as TbMn03 (TMO) 
and N13V2O8 (NVO) which have more than one phase with incommensurately modulated long- 
range magnetic order. Here we give a phenomenological model which relates the symmetries of 
the magnetic structure as obtained from neutron diffraction to the development and orientation of 
a spontaneous ferroelectric moment induced by the magnetic ordering. This model leads directly 
to the formulation of a microscopic spin-phonon interaction which gives explains the observed 
phenomena. The results are given in terms of gradients of the exchange tensor with respect to 
generalized displacements for the specific example of NVO. It is assumed that these gradients will 
now be the target of first-principles calculations using the LDA or related schemes. 

PACS numbers: 75.25,+z, 75.10.Jm, 75.40.Gb 



I. INTRODUCTION 



Recently studies have focussed on a a family of multiferroics which display a phase tran- 
sition in which there simultaneously develops long-range incommensurate magnetic and 
uniform ferroelectric order. The most detailed studies have been carried out on the system 
Ni 3 V 2 8 (NVO)MM A similar analysis of TbMn0 3 (TMO) has also appeared.-^ (For a 
review of both systems see, see Ref. [3.) To illustrate the phenomenon we show in Fig. ^ 
data for the spontaneous polarization as a function of applied magnetic field. 
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FIG. 1: Adapted from Ref. 3. Spontaneous polarization of NVO versus magnetic field applied 
along the a and c axes. 

This data indicates a strong coupling between the magnetic order parameters and the po- 
larization order parameters. To interpret the effect of this coupling, we reproduce in Fig. 
|2] the phase diagram as function of temperature T and magnetic field H applied along the 
crystal a and c-axes. Here we show the three magnetic phases that occur for T > 2K. The 
HTI and LTI phases are two distinct incommensurate phases. A spontaneous polarization 
appears throughout the LTI phase, but does not appear in the other phases. Thus, applying 
a strong enough H to cross the LTI-AF phase boundary will kill the spontaneous polarization 
for H||c or allow it for H||a. Since this transition is discontinuous, the hysteresis shown in 
Fig. [U is to be expected. 

This phenomenon has been explained 3 on the basis of a phenomenological model which 
invokes a Landau expansion in terms of the order parameters describing the incommensurate 
magnetic order and the order parameter describing the uniform spontaneous polarization. 
Already from this treatment it was clear that a microscopic model would have to involve a 
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FIG. 2: Adapted from Ref. U Phase diagram of NVO for magnetic fields up to 8 T applied 
along the a direction (top) and along the c direction (bottom). P (AF) denotes the paramagnetic 
(antiferromagnetic) phase. For -ff||c there is no phase boundary between the P and AF phases. 
The arrowed vertical lines indicate the paths followed in Fig. ^ 

trilinear interaction Hamiltonian proportional to the product of two spin variables and one 
displacement variable. Furthermore, the symmetry requirements of the phenomenological 
model would naturally be realized by a proper microscopic model. In this paper we summa- 
rize the symmetry analysis and briefly describe the microscopic formulation which underlies 
the symmetry analysis. 

II. SYMMETRY OF THE MAGNETIC PHASES 

We first discuss the crystal structure of NVO. The space group of NVO is Cmca (#62 in 
Ref. Isj). Apart from the primitive translations ai = (a/2)i + (b/2)j, a 2 = (a/2)i — (b/2)j, 
and a 3 = ck), the generators of the space group may be taken to be 2 X a two- fold rotation 
about the X - 3jX1S ; Trig, ct reflection taking x into —x, and m z a glide plane which takes z 
into —z, followed by a translation of (b/2)j + (c/2)k. The primitive unit cell contains two 
formula units of NVO and the magnetism is due to the Ni ions (see Fig. |3J). There are 
two crystallographically inequivalent Ni sites which we will refer to as "spine" and "cross- 
tie" sites. The positions of the spine sites within the unit cell are r x = (a/4, d, c/4), r 2 = 
(a/4, — d, 3c/4), r 3 = (3a/4, — d, 3c/4), and = (3a/ '4, d, c/4). The parameter d « 0.136 
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FIG. 3: Positions in NVO of Ni cross-tie sites (filled circles) and Ni spine sites (indicated by + 
and — signs according to whether their y coordinate is larger than or less than that of the layer). 
The dashed lines outline the primitive unit cell. 

indicates that the a-c planes are buckled kagome-like planes. The cross-tie sites within the 
unit cell are at r 5 = (0, 0, 0) and r 6 = (a/2, 0, c/2). 
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FIG. 4: First and second neighbor exchanges paths (respectively X-A-Y or Y-B-Z and X-A-B-Z) 
in the vanadates^ (left) and cuprates (right), where A and B are oxygen ions and X, Y, and Z are 
magnetic ions. In the vandates the first neighbor exchange path is nearly a 90° path, so that the 
nn exchange is weak, whereas in the cuprates the path is a 180 path, so that here Jij J\ ~ 0.1 (see 
Ref. 0), whereas for the vanadates J2/J1 is much larger A 

The Ni spins form weakly coupled chains (or spines) parallel to the a-axis. The nearest 
neighbor (nn) and next-nearest neighbor (nnn) antiferromagnetic interactions (Ji and J2, 
respectively) along the spine compete with one another because, as shown in Fig. in 
contrast to the situation in the cuprates, in the vanadates the nn interaction is anomalously 
weak. A result of this competition is that the spins order in an incommensurate state whose 
wavevector, q, along the a axis is given by- 



cos(ag/2) = -Ji/(4J 2 ) 



(1) 



As the temperature is lowered (at zero magnetic field), the system first orders in the HTI 
phase. Here the spins develop order mainly along the easy (a) axis so that their magnitude 
varies sinusoidally. At a lower temperature the LTI phase is entered in which there is 
additional ordering of the transverse spin components to more nearly satisfy the constraint 
of fixed spin length. 
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FIG. 5: Two adjacent x — z planes in NVO (buckling not shown). Complex spin amplitudes for 
wavevector q (whose real part gives the spin moment). Here 2 X is the axis of two-fold rotation 
about x and m z represents the glide plane which consists of the mirror plane shown at z = c/4 
together with the displacement 6/2 along y (from one plane to the next). 



The above description of the incommensurate phases is too simplistic. Accordingly we 
now turn to a more detailed discussion which takes proper account of the structure within 
each unit cell. The a-component of the spin order at position R (which is the rth site in 
the unit cell) is expressed in terms of Fourier components at wavevector q as 

S ar (R) = [5 aT (q)e^ R + S aT (q)*e-^ R ]/2 , (2) 

where S aT (q) is a complex-valued Fourier coefficient and S aT (— q) = ^^(q)*. In Fig. 
we show the complex quantity S' Q , T (q)e jq ' R whose real part gives the spin order at position 
R. Assuming only a single Fourier wavevector condenses (see below), we see that the spin 
structure is characterized by the six vector Fourier components of the unit cell. Since 
these quantities are three component vectors and each components is complex-valued, this 
structure is characterized by 36 real-valued parameters. As we shall see, the use of symmetry 
drastically reduces the number of parameters needed to characeterize the magnetic structure. 
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Translational invariance indicates that the Landau expansion of the free energy at 
quadratic order in terms of these Fourier components must assume the form 

F 2 = £F 2 (q), (3) 
q 

where 

^(q) = E C ar;a'r';q5' ar (-q)S' Q / T /(q) . (4) 
aa'\TT' 

The eigenvalues of the quadratic form F 2 (q) are the various inverse susceptibilities. We are 
interested in the stability of the quadratic form (since we assume the ordering transition is 
a continuous one). We may plot the least positive (least stable) eigenvalue as a function of 
wavevector q. 
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FIG. 6: Inverse susceptibility as a function of wavevector for a sequence of temperatures for the 
Ji, J 2 model with J1/J2 = 2.56 (as for NVO). Temperature is in units of J\. 

Such plots for a sequence of temperatures are shown in Fig. where one sees that some 
value of wavevector is selected at which the paramagnetic phase first becomes unstable 
relative to the formation of long-range order. Of course, this value of q is determined by 
the interactions between spins, which are not well known. Here we are interested in what 
can be said on the basis of symmetry. A cardinal principle is that we always reject the 
possibility of accidental degeneracy. That is, there is only a single value of q which becomes 
unstable as the temperature is lowered from the paramagnetic phase. (This phenomenon is 
called wavevector selection.) For a ferromagnet, the selected wavevector is q = and for a 
simple cubic antiferromagnetic q = (7r/a)(l, 1, 1), where a is the lattice constant. For NVO 

6 



the situation is more complicated because along each spine one has competing nn and nnn 
antiferromagnetic interactions which cause the inverse susceptibility to have its minimum at 
some nonspecial value of wavevector. This value is^ q ~ 0.28(27r/a)i. The unit cell contain 
several spins, so that the actual incommensurate spin configuration can realize different 
symmetries within the unit cell. In the HTI phase one has incommensurate spin ordering 
mainly on the spines in which the spins are oriented along the easy (a) axis with sinusoidally 
varying amplitude. In the LTI phase, in addition to the order preexisting in the HTI 
phase, significant spin order transverse to the a axis appears, mainly on the cross-tie sites. 
Superficially, this type of structure may seem difficult to characterize. However, by focussing 
on the symmetry of the structure we are able to obtain a convincing phenomenological 
analysis of the symmetry of the magnetoelectric interaction. 

We therefore now turn to the an analysis of the magnetic symmetry. Since the wavevector 
lies along the x-axis, the free energy, F 2 (q) must be invariant under the operations of the 
crystal which leave the x-axis invariant. The group of these symmetry operations includes 
the identity operation S, a two fold rotation about the x-axis, 2 X , the glide reflection through 
a z-plane, m z and 2 x m z . (Some of these are illustrated in Fig. E3) Without going through the 
details of group theory, one can say (in this case) that the eigenvector of the quadratic form 
i*2(q), must also be an eigenvector of the operation 2 X (with eigenvalue either +1 or —1) and 
also of m z (with eigenvalue either +1 and —1). Thus the structure determination is much 
simplified. Instead of having to determine the 6 complex vectors S n of Fig. we use the fact 
that these spin components form a vector which is simultaneously an eigenvector of F%(c\), 
2 X , and m z . Accordingly we may specify (or guess) the four possible sets of eigenvalues of 2 X 
and m z . This set of eigenvalues is technically known as the irreducible representation T, or 
irrep, for short. The actual active irrep is selected as the one which best fits the diffraction 
data. For the HTI phase the active irrep is found to be 2 T 4 for which the eigenvalues of 
2 X and m z are — and +, respectively. For this irrep we have to specify the 3 complex spin 
amplitudes of a single spine site and the y and z components of a single cross-tie site. The 
other amplitudes are then determined by application of the symmetry operators 2 X and m z . 
Since we reject accidental degeneracy, the ordering of the HTI phase can only involve a 
single irrep. In Fig. 

|7J we show the spin function for the HTI phase of NVO corresponding to the eigenvalues of 2 X 
and m z respectively —1 and +1.- Notice that the characterization of the magnetic structure 
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FIG. 7: As Fig. OD Fourier coefficients S aT (q) in terms of symmetry-adapted parameters for the 
HTI irrep T = T4, for which the eigenvalues of 2 X and m z are —1 and +1, respectively. 



is reduced to the specification of only five complex-valued symmetry adapted coordinates. 

When the LTI phase is entered, the situation is much the same. In addition to a spin 
function with the HTI symmetry, we now introduce additional ordering which has a different 
symmetry. This additional symmetry is found to correspond to the eigenvalues of 2 X and m z 
respectively +1 and +1 (of irrep Ti) and the corresponding spin functions (which require 
the specification of four additional symmetry adapted coordinates) are shown in Fig. |H1 

The cystal NVO is inversion symmetric and we have not yet discussed how this affects 
the analysis. Under spatial inversion I the components of spin transform into each other. 
Since spin is a pseudovector, inversion does not change the orientation of the spin, but takes 
it into the spatially inverted position and therefore takes q into —q (which is the same as 
complex conjugation). Also inversion takes the site r x in the unit cell into site r 3 and site r 2 
into site r4. Referring to Figs. [3 and El one see that we have defined the spin components 



rrv^ with appropriate phase factors of i so thal^ 1 * 
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This same relation also holds for the cross-tie order parameters without the necessity of 
introducing phase factors. Now consider the expansion of the free energy in terms of 
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FIG. 8: As Fig. [21 Fourier coefficients S ar (q) in terms of symmetry-adapted parameters for the 
LTI irrep r = T\, for which the eigenvalues of 2 X and m z are both +1. 

these variables, where for notational convenience we set (for each irreducible representa- 
tion) £i(g) = m sx (q), £ 2 (<?) = rn sy (q), £ 3 (q) = m sy (q), and £ 4 (g), and if need be, £ 5 (g) are 
similar cross-tie order parameters. In terms of these variables the quadratic free energy is 



(6) 



n,n' ; r 



where the superscript identifies the irreducible representation. To ensure the reality of ^(q) 
we must have n = c^ n *. Equation (jHJ) implies that 



(7) 



Now for ^(q) to be invariant under X we must have 

JF 2 (q) = £ <*KOTX&(?)] = E = ^(q) ■ (8) 

n,n';r n,n' 

This indicates that c^„/ = c^, n , which together with c^, n (q) = c^ n ,(q)* means that all the 
coefficients <? n n , are real valued. This means that apart from an overall phase factor, all the 
components of the eigenvector are real. So for each representation, T, we write 



£n = ae*r. 



(9) 



where the r n are the components of a real valued unit vector.— Instead of giving the rep- 
resentation T as a superscript on the £'s, 0's, and r's, we will write HTI (for T^) or LTI 
(for Ti). By writing the symmetry adapted coordinates in the form of Eq. (|UJ), one sees 
that this type of incommensurate ordering is characterized by an x-y like order parameter 
<r(q) = crexp(i0) which has a magnitude a and a phase 0. The magnitude of the order 
parameters a are fixed by terms (we have not considered) which are fourth order in the 
spin components. These and higher order terms have no effect on the normalized eigenvec- 
tor as long as one is close to the ordering transition. From Eq. ((7J) we see that for each 
representation 

la = a , 1(f) = -0 . (10) 

Apart from the overall phase, the spin structure of the HTI phase is specified by the wavevec- 
tor q and the 5 real- valued parameters r^ TI and the additional ordering appearing in the LTI 
phase requires specifying 4 real- valued parameters r^ TI and the relative phase 0lti — 4>hti- 



III. MAGNETOELECTRIC COUPLING: PHENOMENOLOGY 

In Ref. |3| it was proposed that one could understand the magnetically induced ferroelec- 
tricity in terms of the following Landau expansion 

F = F M + -^-P 2 + V int , (11) 

where Fm is the free energy of the system in the absence of a nonzero polarization (and 
which is given up to quadratic order by F 2 (q)), P is the uniform electric polarization, and 
Xe is the electric susceptibility, which we assume to be finite, since we assume that when the 
magnetism is absent the spontaneous polarization is zero. Here Vj n t is the magnetoelectric 
interaction which is posited to be of the form 

Vint = Yl C a r;/3r'(q)5 , ar(-q)5'/3r'(q) J P7 , (12) 

where 7 labels the Cartesian components of P. This interaction may be viewed as that 
of an effective electric field proportional to a quadratic combination of spin variables and 
which thereby induces a nonzero polarization only when the magnetic order parameters 
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are nonzero. Because the spin components are expressible in terms of the x-y-\ike order 
parameters introduced in Eq. ©, we may write this as 

V^t = CAB 7 (q)<T A (g)er B (-g)P 7 , (13) 

AB-y 

where A and B are summed over the values HTI and LTI. In the HTI phase the LTI order 
parameter is zero, so that 

Vint = C HT I,HTI )1 {q)\<?HTl{q)\ 2 P 1 ■ (14) 

7 

But since this interaction has to be inversion invariant, the coefficient CHiT,HTi,-y(q) must 
vanish. Thus symmetry does not allow a spontaneous polarization to be induced by magnetic 
ordering in the HTI phase. A simple (but not rigorous) way to reach this conclusion is to 
observe that the phase of an incommensurate wave will vanish arbitrarily close to some 
lattice site. If one then takes this lattice site as the origin, the magnetic structure will have 
inversion symmetry relative to this new origin and hence a spontaneous polarization is not 
allowed. 

In the LTI phase the situation is different. Although the term involving |o"lti| 2 vanishes 
by the above argument, one has 

Vim = Xl jP 7[ C 7(^) Cr HTl(g)crLTl(g)* + C 7 (g)*<T H Tl('?)*CrLTl(g)] • (15) 

7 

Using Eq. (J1(J|) one sees that for this to be invariant under spatial inversion the coefficient 
Cry(q) has to be pure imaginary, so that V in t has to be of the form 

V m t = X! r 7 a HTlCTLTlSin(0 L Tl - 0HTl)-P 7 , (16) 

7 

where r 7 is real-valued. Thus a requirement that for NVO incommensurate magnetic order 
induce a spontaneous polarization is that the two order parameters <t H ti and cr L ti should 
not be in phase, i. e. 0hti 7^ 0lti- An analysis^ of the fourth order terms in the Landau 
expansion of Fm indicates that indeed these two order parameters do not have the same 
phase. We now consider the effect of the symmetries 2 X and m z . From the symmetries of 
the active irreps, (or from inspection of Figs. [3 and |EJ) one can see that 

2z[c H tiClti] = ~~ [""htiClti] (17) 
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and 



m>z [°"htict lti 



o"HTicr LTI 



(18) 



But the product P 7 cthtiClti m ust be invariant under these two symmetry operations. Thus 
P 1 must change sign under 2 X but be invariant under m z . This implies that P can only have 
a nonzero component along b, as is observed-- 3 

To summarize: the Landau symmetry analysis for NVO explains why ferroelectricity 
can only be induced in the LTI phase which is described by two different symmetry or- 
der parameters. (This is reminiscent of the ahalysis of the symmetry of second harmonic 
generation.—) In addition this symmetry analysis correctly predicts that incommensurate 
magnetic order can only induce a spontaneous polarization along the crystallographic b di- 
rection. A completely similar analysis has been developed^ for TbMnC>3 which exhibits a 
similar magnetically induced ferrolectricity.— 

IV. MAGNETOELECTRIC COUPLING: MICROSCOPICS 

As noted in Ref. b[ the form of the trilinear magnetoelectric coupling of Eq. (|12jl leads 
directly to the construction of a microscopic interaction which must underlie the symmetry 
analysis presented above. Since a polarization can only come from an atomic displacement, 
the microscopic interaction must involve the coupling of two spin operators and one dis- 
placement operator. The generalized exchange tensor interaction between spins at sites % 
and j is 



Here we allow for both symmetric and antisymmetric Dzialoshinskii-Moriya (DM) 
interactions . 14 i 15 To get a displacement operator, we simply expand the exchange tensor 
J up to first order in the atomic displacements. Note that a uniform polarization must 
arise from a zero wavevector optical phonon. Since the energies and wavefunctions of the 
optical phonons are unknown, we will express the results in terms of the symmetry adapted 
generalized displacements (GDs) at zero wavevector. We only need to consider GDs which 
transform like a vector under the symmetry operations of the crystal, because only these 
GDs can carry a dipole moment. These GDs are denoted Q an , where a labels the component 




(19) 
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x, y, or z according to which Q transforms. Some of these GDs are easy to construct: a 
GD in which all crystallographically equivalent atoms e. g. Ni spines, Ni cross-ties, or V 
ions, move along the same crystallographic direction a will generate the Q a ,nS for n < 6, 
because NVO has six crystallographically inequivalent sites. The other relevant GDs are 
given elsewhere,—^ but we show a few of these in Fig. El To identify the symmetry of the 
GDs shown here, it is useful to consider the effect of 2 X , m x , and 2 y , which is a two-fold 
rotation about a y-axis passing through a spine site. For instance, consider the 
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FIG. 9: Adapted from Ref. Ilq . Top: generalized displacements of cross-tie Ni atoms which 
transform like y (left) and like z (right). Plain circles represent displacements out of the page and 
a circle with a cross represents a displacement into the page. Bottom: generalized displacements 
of Ni spine atoms which transform like z (left) and like x (right). 

effect of these operations on Q z $. This GD is odd under 2 X (a two-fold rotation about the 
center of the cell), even under m x , and is odd under 2 y . These evaluations verify that this GD 
does have the same symmetry asaz coordinate. This means that, perhaps surprisingly, this 
mode gives rise to a dipole moment in the ^-direction. One can also see that Q z $ couples to 
the GD in which all the cross-ties move in parallel along the z-axis. Suppose the atom-atom 
potential is repulsive. Then when the spine atoms move as shown in Fig. they force the 
cross-tie at the center of the unit cell to move towards negative y and positive z and the 
cross-tie at the corner of the unit cell to move towards positive y and positive z. Thus Q z $ 
couples to a uniform displacement of the cross-tie along the its symmetry label 

indicates. 
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In this development the spin-phonon coupling involves the various derivatives 



(20) 



The algebra needed to construct the spin-phonon Hamiltonian 7^ s _ when the spin operators 
are replaced by their thermal averages is straightforward but tedious and is given elsewhere 
when nn spine-spine, and spine-cross-tie, as well as second neighbor spine-spine interac- 
tions are taken into account.-^* Only terms which involve Q a , n with a = y are nonzero, 
in agreement with the macroscopic symmetry analysis.— To illustrate the results of these 
calculations, we give here the result for nn spine-spine interactions^*^ 



H s - P = 16N uc J2Q yp £ 9f [mgj (g)*mg(g) 

p n,v=x,y,z 

where N uc is the number of unit cells, 3 indicates the imaginary part, and 



(21) 
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- T y pr 
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(22) 



where c = cos(ga/2), s = sin(ga/2), and J a p is the symmetric part of the exchange tensor, 
and D is the Dzialoshinskii-Moriya vector. Note that this result explicitly requires two 
irreps, also in accord with the macroscopic symmetry analysis. It should also be noted that 
due to the low site symmetry, the spin-phonon coupling involves the derivatives of almost 
all elements of the generalized exchange tensor. The next step (which is ongoing) is to 
evaluate the necessary gradients of the exchange tensor from a first principle calculation 
and to incorporate the mode structure and energy of the optical modes at zero wavevector. 



V. SUMMARY 



In this paper I have reviewed the symmetry analysis of the interaction responsible of fer- 
roelectricity in M3V2O8 induced by incommensurate magnetism^ and have also summarized 
some results of the corresponding microscopic theory.^^ The symmetry has also recently 
been applied to some of the phases of ThMC^.^ It seems likely that these approaches will 
prove useful for a number of other mutliferroics. 
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